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Abstract: In composite Higgs models the pseudo-Nambu-Goldstone Boson (pNGB) na¬ 
ture of the Higgs field is an interesting alternative for explaining the smallness of the 
electroweak scale with respect to the beyond the Standard Model scale. In non-minimal 
models additional pNGB states are present and can be a Dark Matter (DM) candidate, 
if there is an approximate symmetry suppressing their decay. Here we assume that the 
low energy effective theory (for scales much below the compositeness scale) corresponds 
to the Standard Model with a pNGB Higgs doublet and a pNGB DM multiplet. We de¬ 
rive general effective DM Lagrangians for several possible DM representations (under the 
SM gauge group), including the singlet, doublet and triplet cases. Within this frame¬ 
work we discuss how the DM observables (relic abundance, direct and indirect detection) 
constrain the dimension-6 operators induced by the strong sector assuming that DM be¬ 
haves as a Weakly Interacting Particle (WIMP) and that the relic abundance is settled 
through the freeze-out mechanism. We also apply our general results to two specific cosets: 
SO{6)/SO{5) and SO{6)/SO{4) x 50(2), which contain a singlet and doublet DM can¬ 
didate, respectively. In particular we show that if compositeness is a solution to the little 
hierarchy problem, representations larger than the triplet are strongly disfavored. Further¬ 
more, we find that composite models can have viable DM candidates with much smaller 
direct detection cross-sections than their non-composite counterparts, making DM detec¬ 
tion much more challenging. 
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1 Introduction 

Despite being an enormous triumph to the Standard Model (SM), the discovery of the 
Higgs boson [1-5] has established the question of how the electroweak scale is stabilized 
under the large corrections from new physics at ultra-violet (UV) scales. The quadratic 
sensitivity of the Higgs mass to UV physics, also known as the hierarchy problem, has 
been one of the leading motivations for searches of new physics at the LHC. One of the 
most investigated solutions to the hierarchy problem is Supersymmetry. Besides providing 
a mechanism for stabilizing the electroweak (EW) scale, the Minimal Supersymmetric 
Standard Model (for reviews see [6-8]) or MSSM also has several attractive features, such 
as gauge coupling unification and viable dark matter candidates. However, current LHC 
searches and the measured Higgs mass impose severe constraints on the MSSM. 

Composite Higgs models [9-14] are also an interesting solution to the hierarchy prob¬ 
lem. Unlike the MSSM, where the Higgs boson is a fundamental scalar, in composite 
models the Higgs doublet is a pseudo-Nambu-Goldstone Boson (pNGB) appearing in the 
low energy theory as a result of the spontaneous breaking of a global symmetry (Q —)■ Ti) 
by a new strong sector dynamics. In analogy to the pions in QGD, the Higgs doublet only 
has derivative couplings and is exactly massless, except for corrections due to the explicit 
breaking of Q. In the Minimal Gomposite Higgs Model (MCHM) [15], based on the coset 
SO{5)/SO{4:), the four pNGBs correspond to the complex Higgs doublet. Nonetheless, 
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other cosets are possible and may contain a higher number of pNGBs degrees-of-freedom. 
It is then an interesting question whether one of the additional pseudo Nambu-Goldstone 
bosons may explain the observed dark matter component of the universe. Since this requires 
at least five pNGBs (the Higgs doublet plus one DM state), one must consider extensions 
of the MGHM. The simplest scenario corresponds to the SO{6)/SO{5) coset [16, 17], which 
provides exactly one additional pNGB. For a model based on the coset SO{7)/G2, see e.g. 
[18]. 

Possible explanations to the DM question can now be tested by an impressive number 
of experiments, possibly bringing us at the edge of discovery. Four complementary ap¬ 
proaches are followed to identify DM as a new WIMP. Direct searches are looking for DM 
scattering off heavy nuclei in underground detectors. For the mass ranges that we are inter¬ 
ested in, LUX is for the time being the most sensitive experiment [19], reaching scattering 
cross-sections of the order of 10“^^cm^ (for masses around 30 GeV). In this field, the next 
data releases are expected for 2015 (LUX-1 year data) increasing the sensitivity by a fac¬ 
tor 5. From 2018, the next generation of experiments with multi-tonne Xenon (XenonlT, 
LUX-Zeplin) should reach cross-sections down to 10“^^cm^. Indirect searches aim at de¬ 
tecting the products of DM annihilation or decay in the form of gamma-rays, neutrinos or 
antimatter cosmic rays. The Fermi-LAT telescope, which is studying the gamma-ray flux 
from multiple astrophysical sources, has provided the first strong exclusion limits on DM 
annihilation cross-sections excluding ((Tann-i^) ~ 3 x 10“^® cm^/s for candidates with masses 
up to ~ 100 GeV for a 100% annihilation into bb [20, 21]. It is very likely that new data 
will soon be available with e.g. AMS-02 for antimatter cosmic rays, Ice-Gube for neutrinos, 
HESS-II for gamma-rays and with a new generation of experiments, such as Gamma-400 
and the GTA. From 2018 onwards, see e.g. [22] for a review, the new generation of indirect 
detection experiments will improve considerably the current results. Collider searches can 
also constrain DM scenarios through missing energy searches, although in most cases the 
constraints are strongly model dependent. In addition, if dark matter has non-negligible 
couplings to the Higgs field, the current measurement of the Higgs decays can also be 
relevant for testing DM scenarios [23-25]. Finally, astro-cosmo probes can provide comple¬ 
mentary constraints, testing e.g. distortions of the Gosmic Microwave Background (GMB) 
signal due to energy injection from DM annihilation in the early universe, see e.g. [26-29]. 

All the data provided by the complementary searches listed above strongly constrains 
dark matter models. In particular, the composite Higgs model 50(6)/50(5) can only 
satisfy all the current DM and LHC constraints if the compositeness scale {F) is > 1 TeV 
and if tudm ^ 200 GeV, as discussed in detail in refs. [16, 17]. Due to the large amount of 
data available from direct, indirect and LHC searches, it is possible to formulate even more 
general statements on the possible composite scenarios, if these are required to provide a 
DM candidate behaving as a WIMP and accounting for all the DM abundance. This is 
the main purpose of this work. Here we consider the class of composite models where the 
only composite states present in the low energy effective theory are the Higgs doublet and 
the DM multiplet, and both are pNGBs. The possible strong resonances and remaining 
pNGBs are assumed to be decoupled in the effective Lagrangian, either because they are 
significantly heavier than the Higgs and the DM multiplet or because they correspond to 
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gauge degrees of freedom, due to a partly gauged Q. We discuss this possibility with a 
minimal knowledge of the UV completion. In particular, we do not specify the fermion 
representations under Q, but parametrize their effects instead. Within this framework, 
any composite Higgs model containing a DM candidate can be mapped at low energies 
(<C A ~ 4ttF) to an effective Lagrangian containing the Standard Model (including the 
Higgs doublet) and the DM multiplet. 

In section 2 we define the class of models we will consider and present the basic 
framework used to compute our results. In section 3 we classify the different scenarios 
according to the DM representation under 50(4) and discuss each case separately. The 
simplest case, consisting of a singlet DM, is discussed in detail in section 3.1. A minimal 
realization of this scenario can be obtained for Q/Ti = 50(6)/50(5), as shown in refs. [16, 
17]. The next simplest scenario consistent with 50(4) C H is the complex doublet DM 
case, discussed in section 3.2. Although we once again present our results in a model 
independent way, we also discuss a realization of the doublet DM case, corresponding to 
the 50(6)/50(4) x 50(2) coset. In section 3.3 we present the results for the triplet and 
general constraints on higher representations. Finally, in section 4 we summarize our results 
and present the conclusions. In the appendix A we define our notation and discuss in detail 
the derivation of the effective Lagrangians used in our calculations. 

2 General Composite Dark Matter 

In composite Higgs models [30-36], the strong interactions responsible for composite¬ 
ness are assumed to break a global symmetry group G down to a smaller symmetry group 
at the F scale. The low energy spectrum, at energies <C A ~ dvrT, consists of the 
massless Nambu-Goldstone Boson (NGB) modes with an effective Lagrangian completely 
fixed by the symmetry breaking pattern Q ^ Ti. For a general G/Ti breaking, the effective 
interactions of the Nambu-Goldstone modes to lowest order in 0{1/F) are given by the 
Callan, Goleman, Wess, Zumino (GGWZ) Lagrangian [37, 38]. In order to avoid large 
corrections to EW precision observables, we assume 50(4) C 77, which implies that the 
CCWZ Lagrangian considered here is 50(4) invariant. 

If one of the NGBs must play the role of the Higgs doublet, G must be explicitly 
broken, resulting in a non-flat potential for the pseudo-Goldstone fields. Here we assume 
the partial compositeness scenario, where the explicit breaking is restricted to Yukawa 
and gauge interactions between the SM and the strong sector. Once these interactions 
are determined, the effective scalar potential for the pNGBs can be computed up to form 
factors describing the strong dynamics. Since the effective scalar potential is generated 
only at loop level, the pNGBs are expected to have masses < F. The fine-tuning in the 
Higgs potential is given by the separation between the F and EW scales. A fine-tuning 
not worse than 1% implies: 

F<3TeV, (2.1) 

but the upper bound in realistic models is usually stronger [39]. Gurrently the strongest 
lower bound on F comes from EW precision tests and direct searches for the strong reso¬ 
nances at the LHG, typically resulting in E > 800 GeV [40]. Although one may argue that 
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incalculable UV effects may hinder the impact of EW precision data on compositeness, at 
this point Higgs data alone can set F > 700 GeV [41] . 

In the MCHM [15], where Q/V. = 50(5)/iS'0(4), the four pNGBs correspond to the 
degrees-of-freedom of the Higgs doublet. Larger cosets can provide a stable DM candidate, 
as shown in refs. [16-18]. In such cases, in order for one of the pseudo-Goldstone bosons to 
be a viable DM candidate, we must assume that both the strong sector and the SM respect 
an exact (or approximate) Z 2 symmetry, so the lightest Z 2 -odd pNGB is (quasi-)stable.^ 
Furthermore, since current LHC data points to a SM-like Higgs, we assume the Higgs to 
be the lightest Z 2 -even pNGB, transforming as a bi-doublet under 50(4). On the other 
hand, the DM candidate can in general be a mixture of distinct SO(4)-multiplets after 
electroweak symmetry breaking (EWSB). Here we only consider the cases where the DM 
state (mostly) belongs to a single 50(4) multiplet, so the mixing with other multiplets 
can be ignored.^ Therefore, the low energy spectrum consists of the Higgs doublet and the 
DM multiplet. We assume all possible remaining pNGBs and strong sector resonances to 
be decoupled in the effective Lagrangian. This approximation is exact in minimal cosets, 
where the Higgs and the DM multiplet correspond to all the pNGBs, as in the particular 
examples discussed in section 3. Additionally, note that higher derivative terms in the G/l-L 
chiral Lagrangian, such as the Wess-Zumino-Witten term [31], can potentially break the 
Z 2 symmetry. Since these terms strongly depend on the coset Q!% and the representation 
of the fermions under Q, which we do not specify, they can not be computed in our model 
independent approach. Hence, we make the simplifying assumption that only those that 
preserve the Z 2 symmetry are present, so the stability of DM is maintained. Furthermore, 
there are particular choices of cosets Q/F. and fermion representations, where the Z 2 - 
violating interactions are zero and the ^ —>■ —^ parity is unbroken, see refs. [16, 17] for a 
specific realization. 

Any composite Higgs model satisfying these assumptions, detailed in the appendix A, 
can be mapped to a low energy effective theory valid for scales <C A, as shown schematically 
by figure 1. Below we discuss the general form of the effective operators relevant for 
computing the DM observables and their impact on the DM phenomenology. 

2.1 Aspects of the Effective Dark Matter Lagrangian 

The most general effective Lagrangian consistent with the assumptions discussed in 
the previous section (see appendix A for more details) contains a large number of unknown 
parameters. The latter correspond to the coefficients of all the possible dimension six 
operators involving SM fields as well as those involving the DM multiplet. Fortunately, 
only a small subset of operators is relevant for computing the DM observables. We also 
consider the minimal scenario, where DM is a cold thermal relic and the observed DM 
abundance is generated through the standard freeze-out mechanism. We only consider tree 
level diagrams in our calculations. In addition, processes where DM annihilates into more 

^Notice that the Z 2 symmetry could already be present within the Q global symmetry group, see e.g.[16]. 

^ In models where, after EWSB, there is a significant mixing between the lightest Z 2 -odd multiplets, the 
lightest Z 2 -odd state is usually charged, so there is no viable DM candidate, see e.g. [42] in the non-composite 
case. 
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Figure 1. Schematic representation of the spectrum for the class of models considered here. At low 
energies (<C F) the effective theory we consider only contains the SM (including the composite Higgs 
bi-doublet) and the DM multiplet. The two possibilities for the DM mass are shown: tudm > rrih 
or moM < ruh- 


than two final states, such as^-|-^ —)■ X + Y + Z are usually phase space suppressed 
with respect to ^ ^ X -|- T. Hence we neglect operators with more than four fields. 

Notice though that in some cases, annihilation into three body final states are well known 
to be important [43-47]. In particular, here we will take into account the contribution 
of annihilations into off-shell gauge bosons C t E -|- . Let us also emphasize that 

dimension five and six operators can become relevant if they contain one or more Higgs 
fields, since after EWSB these operators generate 4-field operators (with H —)• (H) = v). 
Furthermore, the Z 2 symmetry, under which ^ is odd and H is even, requires all operators 
to contain an even number of .^’s. 

Given all the assumptions above, the operators containing more than two powers of ^ 
do not affect the relic abundance or the detection rates. We just need to consider operators 
with zero or two DM fields. Hence, for convenience, we split the effective DM Lagrangian 
in a ^-independent sector and a part containing two powers of ^ 

£eff = /:(°^(^,...)+T(")(F,^,...), (2.2) 

with given by 

£(0)(77,...) =/:sM + /:6, (2.3) 

where Fsm is the SM Lagrangian and Fg contains the relevant dimension-6 operators 
induced by the strong dynamics. Since Fg is independent of its operators are common 
for all the cases to be discussed later. The SM dimension-6 operators have been explored 
at length in the literature [48-53] and are briefly discussed in appendix A. Below we simply 
list the higher dimensional operators considered in our analyses. We use 

£g = ^ - ^\H\^ [{ytQLHHn + ybQLHbn) + h.c.] , (2.4) 
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where a 2 H-, C 4 and Xhg are generated after we integrate out the strong sector, while Ai is 
the dim-4 Higgs self-coupling. As we assume a CP-even Higgs, C 4 is real. The coefficients 
a 2 H,C 4 and Xhg are considered as 0(1) parameters. Notice that the operators in eq. (2.4) 
are weakly constrained by the LHC data [49, 54]. 

In order to reduce the number of parameters that would potentially affect the DM 
observables and simplify our analysis, we consider one single coefficient C 4 multiplying the 
dimension -6 operators involving a direct coupling between the Higgs particle and the top 
and bottom quarks. As far as relic abundance and indirect dark matter detection searches 
are concerned, our choice could be justified by the fact that this extra coupling to top 
and bottom quarks are relevant in different regions of the parameter space. The impact of 
similar operators involving light quarks would be suppressed by the Yukawas (y^). Let us 
emphasize though that this is however not true for direct dark matter searches in which case 
the relevant Higgs-nucleon coupling is dominated by heavy quarks (b, t and c) contributions. 
The latter are due to heavy quark currents coupling to gluons through triangle diagrams 
involving heavy quark loops, see e.g. [55] and references therein^. 

As it is well known, after EWSB the operator proportional to a 2 H generates a non- 
canonical kinetic term for the Higgs field. The rescaling to a physical Higgs introduces the 
factor _ 


R=l/ 


l + 2a2H^ 


(2.5) 


in all Higgs interactions. Even though we include R in our numerical calculations, for 
simplicity we will typically neglect it in the analytic expressions. Indeed i? is ~ 1 for F > 1 
TeV. 

Notice that we do not consider the dimension -6 operators that modify the gauge cou¬ 
plings to fermions, such as (H^a^ D^H) (QlO'^^^Ql), and those which contribute to the 
electroweak parameters at tree level, like H) (W^uB^'^) [48-52]. As it is well known, 

these operators are strongly constrained by the electroweak fit [40]. 

Since we assume that is a S'0(4) ~ SU(2)l x SU{2)ji multiplet, the possible repre¬ 
sentations for ^ are 


(2jL + l,2jn + l) = (1,1), (2,2), (3,1), (1,3), (3,3), (2,4),... . 


( 2 . 6 ) 


In the following we limit our results to the cases 

• Singlet DM: ^ ~ (1,1) 

• Doublet DM: ^ ~ (2, 2) 

• Real Multiplet DM ((^‘^ = 0- ? ~ (^) 1 )j where n = 3, 5,... . 

Our restriction to real representations (except in the doublet case) is justified by the fact 
that the phenomenology associated to complex multiplets would be very similar to the real 

®Let us mention though that the contribution to the amplitude of DM-nudeon scattering related to the 
coupling C 4 is always expected to be a subdominant correction compared to other operators as it appears 
in the following combination: /F'^) (where \(_n is the Higgs-DM coupling). We expect thus 

the impact of choosing a common dimension -6 coefficient C 4 to all quarks to be marginal. 
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multiplet case; the complex multiplet being equivalent to two degenerate real multiplets. 
A real DM multiplet also implies that ^ is not charged under any U{1) symmetry, since 
= C^* = Consequently, ^ is not charged under U{1)y and since the hypercharge 
operator is a linear combination of C SU{2)ii and U{l)x^-, it means that ^ is a singlet 
under SU{2)r. Hence ^ ~ (n, 1). 

Finally we briefly comment on the validity region of the effective Lagrangian used in 
our analysis. The energy involved in the relevant DM processes is typically y/s < 2mj:)M, 
where the upper limit corresponds to non-relativistic annihilation at zero velocity today. In 
the early universe, energies involved are usually slightly larger given that v ~ 0.2 while for 
direct searches the energies involved are much smaller. Therefore, we expect the effective 
Lagrangian approach used here to be valid as long as the new physics states (which could 
affect the DM observables) are well above 2mY>M- In composite Higgs models one naively 
expects the strong resonances to have masses of order A ~ AttF. However, the composite 
partners of the top quark must be lighter in order to avoid re-introducing a large fine- 
tuning in the Higgs potential [39]. Here we assume that the spinorial resonances have 
masses > F, hence our effective Lagrangian is valid for rriDM or <C F, 

which we use as a strict upper limit in our study. It is also important to point out that, 
for rriDM < T, all the operators used in our analysis are far away from any unitarity 
bounds [57, 58]. 


2.2 The Effect of Higher Dimensional Operators 

The specific form of the dimension-4 and 6 relevant operators^ consistent with C 

7i are derived in detail in the appendix A and will be discussed in section 3. It is instructive, 
however, to first comment on the general properties of these operators and how they can 
affect the DM observables. Since we are only interested in 4-field operators after EWSB, 
the relevant dimension-6 operators generated by integrating out the strong sector heavy 
resonances are of the following types 






g'^V„h 




and (2.7) 


where h is the Higgs boson field, represents a SU{2)l x U{1)y gauge field, V' is a SM 
fermion, g is the corresponding gauge coupling, A represents the corresponding dimension-4 
coupling and refers to the typical momentum scale for the process. As all the above 
operators are presumed to have coefficients < 1, one would expect such operators to be 
irrelevant with respect to the corresponding dimension-4 operators for F ^ v. However 
it is often the case that, in order to generate the correct relic abundance, the dimension- 
4 operators coefficients need to be suppressed. In such cases, the dimension-6 operators 
may still play a significant role. Furthermore, the derivative operators can be enhanced if 

~ ^DM > 

^As it is well known in composite Higgs models, it is often necessary to enlarge the Q group with an 
additional U{l)x symmetry in order to obtain the correct hypercharges for the SM fermions [48, 56]. 

®We always mean by dimension-6 operators, the operators suppressed by 1/F^. 
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In order to illustrate the above arguments more concretely, we take the singlet DM 
case as an example. In the non-composite limit (F —?• oo), there is only one relevant 
dimension-4 operator, corresponding to the Higgs-DM coupling (A^/j), also known as the 
Higgs portal coupling. As a result, the parameter space in the non-composite singlet case 
consists only of and rriDM- As it is well known, requiring the DM relic energy density 
— 0.12 [26, 59] fixes the Higgs-DM coupling for a given DM mass and the solutions 
with the correct relic abundance correspond to a line in the X^h-i^dm plane. This is shown 
by the black solid line in figure 2, where we plot the value of the Higgs-DM coupling in the 
non-composite case required to generate the correct relic abundance. As we can see, 

X^i^ does indeed take values <C 1 in most of the parameter space, which may render the 
operators in eq. (2.7) relevant in these regions even though they are suppressed by 1/T^. 
In order to estimate for which values of tudm the inclusion of the dimension-6 operators 
may affect the DM observables, we compare their potential effective strength with respect 
to Xfc. 

Assuming that the coefficients multiplying the operators in eq. (2.7) are 0(1), the blue 
shaded areas in figure 2, for F = 0.8 TeV (left) and F = 2.5 TeV (right), correspond to 
the windows of the parameter space in which these operators can potentially drive the 
relic abundance. In the regions below the continuous horizontal blue line X^f^ < v'^fF'^, 
i.e. the non-derivative operators can be dominant, while the regions below the blue dotted 
curves X^j^ < i.e. the derivative operators can become dominant. As we 

can see, for F = 0.8 TeV, the dimension-6 operators can potentially affect the DM cross- 
sections for niDM ^ 35 GeV. On the other hand, for higher values of F this contribution 
becomes limited to the Higgs resonance, where A|))^ <C 1, or the heavy DM region, where 
the derivative operators can be enhanced. Furthermore, for F = 2.5 TeV, we see that 
the non-derivative operators play no role, except around the resonance. We stress that 
the results shown in figure 2 are purely schematic and are intended only to illustrate the 
potential impact of including the higher dimensional operators in our analyses. However, 
as it will be shown in the next section, the overall features described above do appear in 
our numerical computations. 

It is also interesting to point out that the dimension-6 operators can have very distinct 
contributions to the annihilation and nucleon scattering rates. On the one hand, while the 
momentum transferred in DM-nucleon scattering is C)(keV-MeV), the momentum scale 
involved in DM annihilations in the early universe is of order tudm- Hence, the derivative 
operators which scale as p^ jF‘^ only contribute significantly to DM annihilations and not 
to the direct detection rate. On the other hand, dimension-6 operators involving Yukawa 
type of interactions such as 

[VqQLH^qR) , ( 2 . 8 ) 

where yq is the Yukawa coupling of the g-quark and Cq a coefficient of 0(1), give important 
contributions to the dark matter-nucleon scattering cross section, aoMp, for q = c,b and 
t (through triangle diagrams processes that are, by default, taken into account using dark 
matter tools such as micrOMEGAS, see [60]). 

Let us mention that, in the following, we have considered interactions of the type of 
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moM (GeV) 


Figure 2. Singlet DM model: The black continuous line denotes the values of the DM-Higgs 
coupling in the non-composite case as a function of moM giving rise to the right relic 

abundance. The blue shaded areas correspond to the windows of the parameter space in which the 
dimension-6 operators of eq. (2.7) can potentially affect the DM observables for F = 0.8 TeV (left) 
and F = 2.5 TeV (right). In the regions below the continuous horizontal blue line jF'^ 

while in the regions below the blue dotted curves X^^ < F'^. 

(2.8) with couplings to the 3rd generation of quarks only {Cu,d,c,s = 0)- This assumption 
does not affect the relic density and the indirect detection results, since the contributions 
from light quarks is always suppressed by their Yukawas. The couplings to light quarks can 
however impact aoMp-, especially in the regions where Ag/i is small. Given the values of the 
parameters driving the nucleon matrix element considered in our analysis®, we expect that 
including a coupling to all quarks would at most enhance aoMp by a factor of ~ 4 and only 
if|Ag,|«m2 /F2|C,| (see also [17] for a similar discussion). In order to reduce the number 
of free parameters in our analysis, we have however neglected this possibility. Finally, we 
make one additional simplifying assumption: Ct = Cf,. Once again, this assumption does 
not affect the relic density or direct detection rate, since annihilations to tops or bottoms 
are relevant at distinct regions of parameter space. While (Tdmp could change for non- 
universal couplings, we do not expect that our results are significantly affected by this 
assumption, since we expect Ct Cb<0{l). 

3 Results 

Using the effective Lagrangians derived in appendix A, we now compute the DM ob¬ 
servables and discuss how this constrains the effective operators for a given DM repre¬ 
sentation. The DM relic abundance and other observables have been computed using 
FeynRules [61] and micrOMEGAS [62]. The calculation of the relic abundance includes all 
possible annihilation channels into two body as well as annihilation into off-shell V = W 
or Z boson (^ -b ^ U* -b U), which can be relevant for tudm ^ 'mv- Even though radia- 

®We took the default values of micrOMEGAS3.3 corresponding to = 0.0191, fu = 0.0153 and /f = 
0.0447 in the notations of ref. [60]. 
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tive corrections are already known to affect DM observables in similar non-composite DM 
scenarios, especially the dark matter-nucleon scattering cross section [63-65], it is beyond 
the scope of this paper to evaluate all loop contributions within the composite framework 
(the exception being the well known contributions from triangle diagrams of heavy quark 
loops included by default in micrOMEGAS for the computation of the DM-nucleon scattering 
cross section). One should keep in mind though that such corrections are expected to give 
rise to the lower bound anMp ^ 10“^^ cm^ for non-composite models in e.g. the doublet 
case [64], 

In order to be as model independent as possible, we parametrize the effective La- 
grangian according to the results in appendix A and scan over the 0(1) coefficients of 
the dimension-6 operators in the range [—1,1]. The unknown parameters for effective 
dimension-4 operators (denoted with A couplings) in the DM Lagrangian are also allowed 
to vary freely within a range restricted to [—dvr,dvr]. We also require all the solutions to 
satisfy: 0.0941 < < 0.127. ^ 

In addition to the above conditions, directly taken into account in the scans, we will 
superimpose several constraints from LHC, direct and indirect DM searches. Furthermore, 
the low DM mass region is limited by the Higgs invisible width, which has been recently 
constrained by ATLAS [24] and CMS [25]. Here we use the 95% C.L. upper limit from 
CMS: BR(/i —)• invisible) < 0.58. For the direct detection constraints, we consider the 
latest LUX results [19] on the spin-independent DM-nucleon elastic cross-section as well 
as the projected sensitivity of XenonlT, both at 95% CL. The DM annihilation into SM 
particles in our galaxy could copiously produce stable particles, such as photons, neutrinos, 
positrons or antiprotons. In that framework, indirect detection searches constraints from 
Fermi-LAT [20, 21] and PAMELA [67, 68] can be relevant here. Indeed, a previous analysis 
of composite DM models [17] showed that the bound on the DM annihilation cross-section 
derived from PAMELA measurements of the antiproton spectrum could provide constraints 
complementary to LHC and direct detection searches. The bounds associated to charged 
cosmic rays are, however, well known to suffer from astrophysical uncertainties [69, 70] 
affecting the propagation model. Here we use the Su upper bound constraints derived in 
ref. [70] (figure 6c) for the DM annihilation cross-section into bb. ® Let us mention that such 
a bound has been derived assuming an uncertainty of 50% in the solar modulation of the 
antiproton flux and we take the limits associated to the so called CON propagation model 
(green curve in figure 6c of ref. [70]), which usually corresponds to conservative bounds. 
Let us emphasize that even though similar bounds were derived for other annihilation 
channels such as e.g. —>■ W~^W~, we checked that they provide weaker constraints than 
the direct detection searches. We have thus neglected them. Let us also mention that the 
current preliminary constraints on gamma-ray flux from the Fermi-LAT [21] measurements 
provide even stronger bounds on crv(^^ —)■ bb)\y=o. We have checked that the impact of 
these preliminary results are negligible, so we only consider the constraints from PAMELA. 

^ We used the 2a WMAP 5-year range [66]. For a reference, the Sa WMAP 9-year range is 0.1003 < 
^DMh^ < 0.1273 [59] and the 3 ct PLANCK range is 0.1103 < floAih'^ < 0.1289 [26]. 

®The PAMELA constraints are thus imposed on av{^^ —>■ bb)\y-o rescaling the constraints of ref. [70] by 
the associated branching ratio. 
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The LHC constraints from monojet and monophoton searches are still insufficient 
to compete with the direct detection or invisible Higgs decay searches. Nonetheless, for 
both the doublet and triplet representations, the DM multiplet contains charged states, 
which can decay to DM and be detectable at the LHC. The main final states are then 
missing Et plus VL’s, Z's and/or Higgses. Using SModelS [71] we were able to recast 
the LHC constraints for SUSY searches and apply them to the scenarios considered. We 
have explicitly computed the relevant production cross-sections and verified that the LHC 
constraints do not lead to any new excluded region. Therefore we will not consider them 
in the following results. 


3.1 Singlet DM 

3.1.1 The Generic composite singlet DM 

The singlet DM case corresponds to the following effective Lagrangian (for details see 
appendix A.l) 


iff?" - L (i + 


2 ^ 2 
{vtQLHHR + ybQLHbji) + h.c. 




(3.1) 


i?2 


resulting in the following DM mass after EWSB: 


+ -^3 j y • 


Aq A 


(3.2) 


The coefficients adi,d 4 ,X'^ are taken to be real 0(1} parameters, while the A 3 coupling is 
allowed to vary in the window [—dvr, dvr] and /rg is the DM bare mass. Except for a^i, which 
is fixed hy Q/Ti, all the other coefficients depend on the fermion embedding in Q. Therefore 
in our model independent approach we take these to be free parameters. As discussed in 
Sec. 2 . 2 , we neglect dimension -6 operators involving light quarks and we assume a universal 
coefficient (^ 4 ) for the top and bottom couplings. 

Eor ruDM < rnt, the DM observables are mostly determined by the effective ^(pi) — 
^(P 2 ) — h(ph) coupling, which in the composite case is given by: 

V = Q-“^i|l). (3-3) 

where 

A = As ^1-b Agl^^ (3.d) 

and pf^ is the (off-shell) Higgs momentum. 

It is important to point out that the effective DM-Higgs coupling relevant for the direct 
detection (DD) rate and the relic density (RD) (or ((Jann-'w)) can take very different values, 
due to the momentum dependence in eq. (3.3). While for annihilations through the Higgs 
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s-channel we have ~ for direct detection we have t-channel scattering, hence 

~ 0. Therefore it is convenient to define the effective couplings in each of these regimes: 


^ih\ph^2mr>M — ~ r, 


A 


m 


DM 


2 

= Ag/j = — , 


F 2 


(3.5) 


where A|^ is the relevant coupling for the calculation of the relic abundance (DD 

cross-section®). 

In the low mass region, where annihilations to tops are closed, both the DM-nucleon 
spin independent scattering cross-section, aoMp, and the thermally averaged annihilation 
cross-section times the relative velocity at the freeze-out time, (u ann .u). are proportional 
to the Higgs-DM effective coupling: 


CTDMp OC ^ 


m 


DM 


i ?2 


A" ^ 


(3.6) 


In analogy to the non-composite case, once we impose Q.DMh'^ — 0.12, is hxed for 
each value of mnM- Denoting by A^*^ the value of the coupling required to generate the 
correct relic abundance in the non-composite case (shown by the black curve in figure 2 ), 
we have: 

= (1)'- (Ar±4«=|:£)' 

/ 777,^ \ ^ 

=>(yDMp OC ( xfh ± 4arfi j or (tdmp/(^dmp 

where is the DM-nucleus cross-section in the non-composite case. From the above 

result we see that ctdmp can be either enhanced or suppressed with respect to its non¬ 
composite value, depending on the values of a^i and itidm- Let us stress that the above 
results are only approximate, since we have neglected the contributions from ^ 4 , both for 
relic abundance and direct detection considerations. The impact of is, however, included 
in our numerical analysis. 


= 1 ±4 


Orfl '^dm 




i?2 


(3.7) 


3.1.2 Scan results 

We now present the full numerical results for the DM observables for a random scan 
over the coefficients of the effective Lagrangian. Since in some cases the correct relic 

®Let us mention that, since we have not considered DM contact interactions with light quarks, the di 
coupling has only a minor impact on the general DD picture. If couplings to light quarks were included, the 
results could be modified in some small regions of parameter space, where |A^rr| ^ |d 4 |m^/F^, see Sec. 2.2 
for a brief discussion. 
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abundance requires extremely small values of A 3 , A 3 and Xh6, we scan logarithmically over 
their allowed range. All the other parameters are scanned linearly within the ranges: 


10 GeV < niDM < F , 

—dvr < A 3 < dvr, 

-1 < A'3 < 1, 

10-6 <\hq<1, 

— 1 < C 4 , a 2 _H', ^ 4 , fldl < 1 ■ (3.8) 



moM (GeV) (GeV) 


Figure 3. Generic composite singlet DM: Values of the DM annihilation cross-section at zero 
velocity for F = 0.8 TeV (left) and F = 2.5 TeV (right) satisfying to 0.941 < < 0.127. The 

channels giving rise to the largest branching ratio to the annihilation cross-section are shown with 
different colors. We also show as a solid gray line the value of (uann-r’) for the non-composite case 
{F —>• 00 ). 

First, in figure 3, we show the results for the DM annihilation cross-section ({a »„n -v)) 
today at zero velocity as a function of mp)M- For comparison purposes we also show the 
respective values for {(T ^„n -v) in the non-composite case (solid gray line). As we can see, 
for both F = 0.8 TeV (left plot) and 2.5 TeV (right plot), the annihilation rate follows 
closely the non-composite case, as expected from the relation oc 1/((T ann .u). We 

notice, however, that while in the F = 2.5 TeV case the dominant annihilation channels 
are the same as in the non-composite case (^^ —)• bb or —>■ W~W~^), for F = 0.8 TeV, 
both annihilations to tt and hh can be dominant at large DM masses. This is mostly due 
to the fact that, for small values of F, the non-derivative operators can be 

comparable to X^h, as schematically shown in the left hand (LH) side plot of figure 2. On 
the other hand, for F = 2.5 TeV, these operators are always suppressed. 

The DM-nucleon scattering cross-section {aoMp) is shown as a function of uidm in 
figure 4 for F = 0.8 TeV (left) and 2.5 TeV (right). All points give the correct DM relic 
abundance. Once again the non-composite case is shown by the solid gray line. Going 
from low to high rrioM it is possible to recognize in figure 4 some known features [16, 17]: 
the Higgs resonance region {tudm = nT./i/2), the —)■ WW threshold {rriDM — raw) and 
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the high mass region {moM ^ rnt). We also show by distinct colors the absolute value of 
the coefficient adi for each point. As we can see, the low mass region {ttidm ^ 50 GeV) 



Figure 4. Generic composite singlet DM: Values of the DM-nucleon scattering cross-section as 
a function of moM- All points satisfy 0.0941 < < 0.127. The absolute values of the 

dimension-6 coefficient adi are shown with the gradient of colors. The LUX 95% CL bound is also 
shown. 

is unaffected by the dimension-6 operators induced by the strong sector. This is already 
expected from the discussion in section 2.2, where we have shown that the low mass region 
is mostly insensitive to these operators. 

Once rriDM ^ 50 GeV, aoMp can take values drastically different from the non¬ 
composite case. As expected from eq. (3.7), this is due to the contribution of the derivative 
operator proportional to adi (at least for tudm ^ ^t)- Indeed, given that we can get the 
right abundance with ~ i-®- for small values of A, ctdmp can become 

very suppressed, as seen in figure 4. In addition, since adi can take both negative and 
positive values, we can also have the correct DM abundance with A larger than in the 
non-composite case, resulting in an enhancement of aoMp, as shown by the points above 
the gray line in figure 4. 

For moM > annihilations to top pairs become kinematically allowed and (cJann-'c) 
receives a contribution from the operator proportional to d^. In this case, we can once 
again generate the correct relic abundance, even when \^h ^ ■ This corresponds to 

all the points at large masses with small values of ajpMp and adi in the LH side plot of 
figure 4. Notice though that this feature is not as strongly present in the right hand (RH) 
side plot of figure 4 for F = 2.5 TeV. This is simply due to the fact that for such large 
F values and moM > fnt, non-derivative operators, such as the one proportional to d^, 
are always subdominant, as already anticipated by the RH side plot of figure 2. Finally, 
we also emphasize that for F = 2.5 TeV, (Tdmp only deviates significantly from its non¬ 
composite value for higher DM masses {moM ^ 150 GeV), as expected from eq. (3.7) and 
in agreement with the discussion in section 2.2. From figure 4 we already see that for 
m-DM > 200 GeV to be consistent with LUX we need la^il <0.2 for F = 0.8 TeV, higher 
values are however possible if we increase the scale of compositeness. 
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3.1.3 Experimental constraints 

We now discuss the current and future constraints on the generic composite singlet 
DM scenarios. As seen in figure 4, the direct detection cross-section can be significantly 
suppressed or enhanced in the composite case. Nonetheless, the low mass region still 
behaves mostly as Higgs portal DM, being strongly constrained by direct, indirect DM 
searches and invisible Higgs decays. 

We summarize the constraining potential of each of the experiments mentioned at the 
beginning of section 3 in the planes A — niDM and — ttidm of figure 5. Let us hrst 
mention that, in some regions of the parameter space, very small values of A are necessary 
to account for the right relic abundance. Although this is typically consistent with most 
experimental constraints, in the composite framework, both Ai (the Higgs self-coupling) 
and A are expected to be generated by the explicit breaking of G and to be of similar size, 
i.e. 0(Ai) ~ 0(|A|) ~ 0.1. An explicit example can be found in Ref. [17], where A ~ 0.13 
for G/Ti = SO{6)/SO{5). Therefore we consider the region with |A| < 10“^ as theoretically 
disfavored. 
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Figure 5. Generic composite singlet DM: Experimental constraints in the |A| — moM (top) and 
Odi — tudm (bottom) planes. We show results for F = 0.8 TeV (left) and 2.5 TeV (right). See text 
for details. 

In figure 5, the top (bottom) plots illustrate, in the plane |A| — moM (^di ~ '^dm)-, 
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for F = 0.8 (left) and 2.5 TeV (right), which regions are currently excluded by direct 
DM searches (with LUX in red), by indirect DM searches (with PAMELA in magenta) 
and due to large contributions to the Higgs invisible width (in yellow). The pink regions 
are currently allowed and will eventually be probed by XenonlT. The gray regions are 
experimentally allowed by current data and outside the planned sensitivity of XenonlT, 
but they are theoretically disfavored. In the bottom plots (for ) most of the gray points 
are hidden behind the region that will be probed by XenonlT. 

As expected from the discussion above, the low mass region {niDM ^ 100 GeV) is 
almost entirely excluded either by indirect, direct searches or the Higgs invisible width 
constraints. Let us emphasize that this is valid for both small and large values of F. The 
only points still allowed in the low mass region lie around the Higgs resonance and will 
be tested by XenonlT. However most of these points have |A| < 0.01, which is difficult 
to realize in specific models without tuning. Interestingly, for ttidm > 200 GeV, a large 
portion of the parameter space (blue points) with |A| >0.1 evade all current experimental 
constraints as well as the projected XenonlT sensitivity, in contrast with the non-composite 
Higgs portal scenarios. 

Let us emphasize though that the allowed values of adi are already highly constrained 
by LUX and have good prospects to be tested by XenonlT. Since in a specific model 
(where G/Ti is specified), o^i is no longer a free parameter, the allowed region of parameter 
space will be strongly constrained by direct DM searches. For example, in the case of 
QjF. = iS'0(6) X SO{5) discussed below, adi = 1/2 and a large portion of the parameter 
space is already excluded for low values of F. For F = 2.5 TeV, if no DM signal is observed 
in XenonlT, there will be most probably only a small range of DM masses allowed for a 
given adi value. 

3.1.4 Specific Case: G/F = SO{6)/SO{5) 

Although we have taken adi as a free parameter, in a specific composite model, adi is 
completely fixed by the choice of the coset G jF. Therefore it is instructive to discuss how 
the general results presented above apply to a specific model. The minimal coset consistent 
with F C SO{F) is G/F = SO{6)/SO{5), where the 5 pNGBs correspond to the Higgs 
bi-doublet and the DM singlet. This case has been discussed in detail in refs. [16, 17], so 
we limit our discussion to re-interpreting the results obtained previously according to the 
general parametrization introduced in eq. (3.1). 

The GIF = SO{&)/SO{f)) case corresponds to the particular choice of values: 

a2H = adi = 1/2 . (3.9) 

Notice that we keep C 4 , ^ 4 , A 3 and A 3 as free parameters in our analysis. Once adi 
and a 2 H are fixed, we can derive from eq. (3.7) a simple expression for the enhancement 
(or suppression) of croMp with respect to its non-composite value: 

aoMp/a^Mp = (1 ± rf , where r = , (3.10) 
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where once again is the value of the Higgs-DM coupling in the non-composite case 
required to produce the correct relic abundance. Unlike the general case discussed above, 
now the enhancement or suppression factor r is fixed (for a given value of moM)- Hence, 
instead of spanning a wide range of values, ctdmp now lies in two possible branches, one 
for the case 1 -|- r (enhancement) and one for 1 — r (suppression). 



Figure 6. Singlet DM with Q/'H = SO[%)/SO{^): Values of the DM-nucleon scattering cross- 
section as a function oimpyM- All points satisfy 0.0941 < VlpyMh^ < 0.127 while the black continuous 
lines correspond to the approximation of eq. (3.10). We also show with different colors the channels 
contributing with the largest branching ratio to the effective annihilation cross-section relevant at 
the time of freeze-out. 

In figure 6, we plot (Jdmp as a function of moM, but now for the particular case of 
QjH. = 50(6)/iS'0(5). We also show with black lines the approximation of eq. (3.10) for 
the enhanced {(JuMp = ^oMp suppressed [aoMp = ^oMp values of 

the DM-nucleus cross-section. As we can see, the solutions obtained through the random 
scan follow closely the approximation of eq. (3.10). It is interesting to notice that for 
F = 0.8 TeV all of the enhanced solutions have already been excluded by LUX, while only 
the low mass region {uidm ^ 120 GeV) of the suppressed branch is allowed. The situation 
is however quite different for higher values of F, where we see that most of the parameter 
space for both branches has not yet been probed. 

Finally we point out that, for F = 0.8 TeV, we find no viable solutions for rrioM > 500 
GeV, even before applying the direct, indirect and Higgs experimental constraints. Once 
i^DM > iT^hi the annihilation into two Higgses can contribute significantly to ((J aTin .u), 
which can be large even if ~ 0 due to derivative coupling between DM and Higgs [16]. 
As a result, the annihilation cross-section becomes too large, suppressing the relic abun¬ 
dance below the WMAP value. However this only happens for small enough values of F, 
as seen in figure 6. 

3.2 Doublet DM 

If we consider non-minimal cosets, such as e.g. G/Ti = S0(6)/S'0(4) x 50(2) (see 
section 3.2.3 and e.g. ref. [56]), it is possible to have non-singlet representations for the DM 
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multiplet. Models with multiplet DM differ greatly from the singlet case in two aspects: 

• they include tree level (unsuppressed) DM-gauge boson couplings; 

• they allow for co-annihilations between the DM multiplet components. 

Since we require invariance under the simplest DM multiplet representation corre¬ 

sponds to a complex doublet. In this case, the non-composite limit {F —)> oo) corresponds 
to the Inert Two Higgs Doublet Model or IDM [72]. 


3.2.1 Generic composite doublet DM 

The generic doublet case Lagrangian corresponds to eq. (2.4) with the addition of 


= (A.?)’ - MekP - As il + - A4 il + ^|/fp) If'Hp 




hY + h.c. 


-f 


Qdl 

2T2 




(3.11) 


+ ^ + h.c.) + h.c.) + ^ [d^ + h.c.)]' 


-f 


OdA 
^ [' 














^ 1 ^ 1 ^ {vtQLHHR + UbQLHbji) + h.c. 

de 

Y2 [' 


{vtQLdHHji - ybQLdHbn) + S^QLaH%R + yt£}ai'^QiaHtR + h.c. 


where the coefficients a(ii,..,,d 4 ) <^ 4 ; ^ 6)-^3 4 5 are taken to be real 0(1} parameters, the cou¬ 
plings A 3 ^ 4^5 are allowed to vary in the window [— 47 r, 47 r] and is the DM bare mass. The 
Pauli matrices are denoted by ct = {cJi, cr 2 , 1 T 3 } and (j)^ D Y ^ “ (-DmV')V- We also 

define 




(3.12) 


4/2 ( y ’ 

see section A.2 for more details. After EWSB, the DM multiplet acquires the following 


masses: 


rnoM = 


= 1^1 + y (a 3 + A4 + As + 
mlo = ^ (a 3 + A4 - As + M3±A4^^^^2 

= /^| + (l + ■ 


(3.13) 

(3.14) 

(3.15) 


In the results presented below we consider only the cases where rujjo < m^o, rn}j ±, so 
is always the DM.^® From eqs. (3.13)-(3.15) we see that the three masses are in principle 
independent. However, large mass splittings are limited by the perturbativity requirements 
on the Aj,A' couplings. 


^"However we expect all the results derived below to remain valid in the case is the DM particle. 
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Although the inclusion of the dimension -6 operators induced by the strong sector 
introduces 9 new coefficients, most of the DM properties depend on a subset of these (or 
appropriate linear combinations). For instance, the Higgs-DM coupling is given by: 

= 2 ~ ^0,d2 + dOcis) ^d3~ - p2 ’ (3.16) 

where ph is again the Higgs momentum and 

A= J;a.(i + A'|1) . (3.17) 

7=3,5 ^ ^ 


It is useful to distinguish between the DM-Higgs coupling which enter s-channel annihila¬ 
tions {p1 ~ 4m|)^) and t-channel scattering (p^ ~ 0). The first case is relevant to (a »r,u -v), 
while the latter enters the calculation of croMp, hence we define: 

^^ h \ pl = 4 nil ^^ = = 2 “ ^“<^1 + ^“<^2 + 20 ^ 3 ) , 

kh\pl=o = = ^ - 2a,3^ . (3.18) 


While is very similar to eq. (3.3) for the singlet case, the coupling relevant for direct 
detection now contains a contribution from The latter, however, is only relevant for 
large DM masses. Also notice that the vertex is not modified by the dimension- 

6 operators driven by the 0^4 or 0^2 coefficients, since the corresponding operators vanish 
for C ^ (0, H°)/V2 and H^{0,h + v)/V2. 

In the large mass region {niDM S> mw) the annihilation to gauge bosons is open 
and become the dominant annihilation channel. Even in the pure gauge limit, i.e. when 
all quartic couplings are set to zero, ((J ann .u) always gets a non-zero contribution from 
—)■ WW, ZZ. Let us focus on the analysis of the latter channel, neglecting the co¬ 
annihilations, and the —?■ hh, tt channels. The annihilation cross-section into gauge 

bosons in the high mass region can be approximated by: 


{cra.nn-v)\HoHo^VV — 

Mt VtVt) ^ 

Ml ^ VlVl) 


32Trml)M 

4(1 + 


(2\Mt\‘^ -I- \Ml\^) , where 
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/ Am? ti2 ^ 

+ 2arf2^ 
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DM 


,(3.19) 


where Mt and Ml are the contributions to the amplitude from transverse and longitudinal 
modes respectively and gv = gig/cw), = {j^dm + ^h+{ao)) = "^H+iAO) ~ 

for V = W{Z). From the above expressions we have: 


(cTann-v) > 


Wirm^M 


[AdT’P — 


9v 


327rm 


2 ’ 

DM 


(3.20) 
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resulting in too small relic abundances, unless rriDM ^ 500 GeV. This is a well known result 
from the IDM [42, 73] and it still holds in the composite case, since the contributions from 
the dimension -6 operators to Mt are negligible. Notice that, in contrast with the non¬ 
composite case [42], for sufficiently large moM the annihilation cross-section is driven by 
the term proportional to 0 ^ 2 - Let us stress though that, as mentioned in section 2, our 
effective Lagrangian approach is only valid for rriDM ^ F- 

From eq. (3.19) we see that if the A4 l contribution had to be suppressed one would 
need small mass splittings among the components of the DM doublet <C 1) 

and the combination 

—2 -(3-21) 

should be small. Since can not be made arbitrarily large (otherwise —)• hh 

enhances (o'ann-r')), we will see below that the above condition limits the possible values of 
ad2- 

3.2.2 Scan results and experimental constraints 

For the numerical results presented below we scan over all the parameters in eq. (3.11) 
plus the ones in the Higgs effective Lagrangian (eq. (2.4)) within the following ranges: 

2 GeV < rriDM < F , 
rnDM < rriAO < F , 
max(70 GeV, m^M) < < F, 

—47r < A < 47r, 

10-6 < Ah6 < 1, 

— 1 < A 3 , A 4 , A 5 < 1, 

— 1 < C 4 , a2H, di, de < 1 ) 

— 1 < «dl, 0(^2, “dS, ad4 < 1 • 

Again, we scan logarithmically over A’s ranges and linearly on all the other parameters 
range. Notice that we have conveniently replaced Ai, A 2 , A 3 and //| by moM, fnu± 

and A. We have also imposed on the scans that viable models have a stable EW vacuum 
with {^) = 0 [74] and are consistent with LEPII constraints on BSM states, following 
ref. [75]. We require that new charged scalar particles {H^) have masses larger than 70 
GeV, see [76-78]. Eurthermore, as it is well known from IDM results [73], for moM ^500 
GeV the correct relic abundance can only be obtained for small mass splittings within the 
DM multiplet. In this mass region we have imposed that mfj± — rriuM < 30 GeV and 
ruj^o — rriDM < 30 GeV.^^ We can divide the discussion into two mass regions: a low mass 
region for tudm < mh and a high mass region for mDM > rnh- 

We have checked that the mass splittings obtained in the scans are always below 35 GeV for all F 
scales assnmed here. For the purpose of this work, the range used in the scans (< 30 GeV) is sufficient to 
discuss the composite DM phenomenology. 
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Figure 7. Generic composite doublet DM: Values of the DM-nucleon scattering cross-section 
(eroMp) as a function of todm satisfying 0.0941 < < 0.127. In the top plots {F = 0.8 

and 2.5 TeV) and the bottom LH plot [F —)■ oo), we show the low mass region. The channels 
contributing to the largest branching ratios to the effective annihilation cross-section at freeze out 
are depicted with different colors. In the bottom RH side plot, we show <tdmp in the high mass 
region for F = 0.8 TeV (dark blue), 2.5 TeV (light blue) and the non-composite case (F —oo in 
gray) superimposed. 


In figure 7, we show the DM-nucleus scattering cross-section for F = 0.8 TeV and 2.5 
TeV for ttidm < W'h in the top plots. For reference, we also show the corresponding non¬ 
composite limit, F —)• oo, in the bottom LH plot. In this region of the parameter space, 
LEPII [75, 76] constraints require a large mass splitting between and or H^. In 
particular, all the region with itidm + nT,^o < mz is excluded due to the constraints on the 
Z-invisible width. This combined with the lower bound on m^± makes co-annihilations 
negligible for tudm ^ 60 GeV. Since annihilations to gauge bosons are kinematically for¬ 
bidden in this region, the main features distinguishing the doublet and singlet cases - 
co-annihilations and DM gauge couplings - are negligible for small DM masses and the 
doublet DM behaves effectively as singlet DM. 

Once rriDM ^ the —)• VV process rapidly enhances the DM annihilation 

cross-section.^^ For mpiM ^ HO GeV this enhancement is still partially suppressed by 

^^Since we include the effects of annihilation to off-shell IT’s and Z’s (H^H^ —> V*V), this channel 
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kinematics and can result in the correct relic abundance if there is destructive interference 
between the —)• h* —)■ VV and —)• VV amplitudes [79], which happens for 

A|^ < 0. For larger values of the DM mass, the contribution from annihilation into gauge 
bosons is no longer suppressed and generates too small relic abundances, independently 
of the choice of parameters. This is a well known feature of the IDM and remains valid 
in the composite scenario. Nonetheless, in the region where —?■ VV is (partially) 

suppressed, the contributions from the dimension-6 operators to A|^ allow for a wider 
range of A values consistent with — 0.12. As a consequence we obtain larger viable 

parameter space in figure 7 for F = 0.8 TeV and F = 2.5 TeV (top plots), when compared 
to the non-composite case (bottom LH plot). Also, the allowed low mass region extends 
to slightly higher values of tudm due to the presence of the dimension-6 operators. 

In the bottom RH side plot of figure 7 we show values of aoMp in the high mass 
region. We see that viable solutions are only obtained for moM ^ 500 GeV. The relative 
mass splittings between the dark components has to be small. Typically we obtain Am^ < 
20 GeV <C rriDM- Since 0^2 can take negative values, the allowed mass splittings in the 
composite case are slightly larger than in the IDM, due to the cancellation between the a(i 2 
and Am| terms in eq. (3.19). In addition, in order to properly quantify the correlations 
between A^ and ad 2 , we show in hgure 8 the values of ad 2 with gradient colors in the 
A|^ — rriDM plane for F = 0.8 and 2.5 TeV. We see that 0^2 and A|^^ typically have 
opposite signs guaranteeing a suppression of the combination (3.21). We thus see that 
all the contributions to A4 l need to be suppressed in order to account for the right relic 
abundance. For F = 0.8 TeV, all points have \ad 2 \ ^ 0.3, while for F = 2.5 TeV \ad 2 \ can 
be as large as 1, due to the suppression factor v'^jF'^ ~ 0.01. 



^DM (GeV) (GeV) 


Figure 8. Generic composite doublet DM in the high mass range: Values for the dimension-6 
coefficient 0^2 in the Xf^-rriDM plane for F = 0.8 TeV (left) and 2.5 TeV (right). 

From hgure 7 we see that aoMp can take larger values in composite scenarios {F = 0.8 
or 2.5 TeV) than in the non-composite case {F —)> 00 ). This is particularly visible in 
the large mass range in the bottom RH side plot where we have aoMp ^ 10“^'^ (10“^^) 

becomes relevant for slightly smaller values of the DM mass {moM ^ 60 GeV), as seen in figure 7. 
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cm^ for F ^ oo {F = 0.8 TeV). This is easily understood looking at eq. (3.18). In the 
non-composite case = X^/f = A/2 and aoMp oc A^, while in composite scenarios we 
have 

+ 2a,2) . (3.22) 

As a result, |A|/^|, the quartic coupling driving aoMp, can get extra contributions due to 
the adi terms in composite models. 

Figure 9 illustrates the experimental constraints on the generic composite doublet 
case, the color code is the same as the one used in section 3.1.3. Once again, points 
with |A| < 0.01 are theoretically disfavored in realistic composite models, where we expect 
O (Ai) ~ O (Xi) ~ 0 . 1 . Similarly, unless the 0^3 and A contributions cancel each other 
in eq. (3.18), at the price of some fine tuning, we will typically have O (A) O (AjT) 
and too low values of X^f^ will also be disfavored. Hence, we consider all points with 
min(|A|(P|, |A|) < 10“^ as theoretically disfavored and they are shown in gray. As we 
can see, in the low mass region all points are either excluded by current experimental 
constraints or involve too small values of |A|^^|, |A| to be theoretically acceptable, except 
for a small number of points (with min(|A^|, |A|) > 10“^) lying between the resonance 
region and the W threshold for large values of F . Therefore, most of the low mass region of 
generic composite doublet DM models is already experimentally excluded or theoretically 
disfavored. In contrast, the high mass region has just been marginally probed by LUX and 
it even XenonlT will not be able to test it in its entirety. 




Figure 9. Generic composite doublet DM: Experimental constraints in the |A|/^| — rnoM plane for 
F = 0.8 TeV (left) and 2.5 TeV (right). See text for details. 


3.2.3 Specific Case: g/U = SO{<o)/SO{F) x 50(2) 

Once again it is interesting to apply the general results presented above to a specific 
model. Here we take g/% = 50(6)/50(4) x 50(2), where the only pNGBs are the Higgs 
bi-doublet and the complex doublet DM. This case corresponds to the effective Lagrangians 
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in eqs. (2.4) and (3.11) with: 


a 2 H = 2 > “<^2 = 1 and a^i = 0^3 = 0^4 = 0. (3.23) 

Notice that, according to the results of figure 9, we do not expect many allowed solutions in 
the low mass range after the LUX constraints have been imposed for |A| > 10“^. Therefore 
here we only discuss the high mass {moM > 500 GeV) solutions. 

One important consequence from specifying the 01% coset is that it hxes the value of 
0 ^ 2 - As we have seen above, the high mass region requires either large values of F or small 
values of 0^2 (see figure 8 ). Hence, for 0/% = S'0(6)/S'0(4) x 50(2), where 0^2 = Ij we do 
not expect viable solutions for large DM masses and small values of F. This is highlighted 
in the LH of figure 10, where we allow to vary outside the 0.09 < ^ 0.12 

range and show the relic density as a function of ttidm for different F values, shown by 
the color gradient. For this plot, we take m^o = = uidm and A = A' = 0 since this 

choice approximately maximizes (minimize) the relic abundance ((cJann-n)) and we take 
= d-Q = 0 , since these coefficients play no major role in the high mass region. 

We can understand the behavior of in the LH of figure 10 by analyzing the 

rriDM and F dependence of (o-ann-'y)(?? —^ VV) derived in eq. (3.19), but in the “pure 
gauge limit” (all quartic couplings set to zero). The relic abundance first grows with 
rriDM due to the usual dependence oc 1/(cT ann .u) oc This is indeed what 

we observe for the lowest range of masses in figure 10 (LH side). The non-composite case 
in the same limit is plotted for reference with gray points and also show this growing 
dependence in moM- When tudm becomes comparable to F, the dimension -6 derivative 
operator starts to dominate (cTann-'w) and the contribution from the 0^2 term results in 
OlDuh'^ oc l/(o' a.nn .t') OC For a given F scale there is thus a hxed mass range, 

bounded from below and above, of viable DM candidates, i.e. in agreement with the CMB 
bound depicted with vertical dashed gray lines in figure 10 (LH side). This viable mass 
range widen with increasing F. Indeed, we see that for e.g. F = 2.4 TeV we have 550 GeV 
< ruDM ^ 1-2 TeV, while for e.g. F = 2.1 TeV we have 600 GeV < moM ^ 800 GeV. In 
addition, once F <2 TeV, Q.DMh'^ < 0.09 for any value of tudm in the large mass range. 
Therefore, for viable composite scenarios with 0/% = 50(6)/50(4) x 50(2) one needs to 
consider F > 2 TeV. 

In figure 10 (RH side) we show results for a scan on the parameter space of 01% = 
50(6)/50(4) X 50(2), where the non-derivative coefficients are still free parameters. We 
have explicitly verihed that even when all non-derivative parameters are allowed to vary 
there are no viable solutions for F = 0.8 TeV. We thus present points only for F = 2.5 
TeV and F —)• 00 . As shown by the results, for F = 2.5 TeV, we only obtain the correct 
relic abundance for 500 GeV < muM ^1-1 TeV, as expected from our discussion above. 


3.3 Triplet DM and Higher Representations 

In sections 3.1 and 3.2 we have discussed the cases of the singlet and complex doublet 
DM, which can be realized for specihc choices of the coset 01%^ such as the 50(6)/50(5) 
and 50(6)/50(4) x 50(2) ones discussed above. If one wants to consider even larger 
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m^M (GeV) 


Figure 10. Doublet DM with QjT-L = SO{<6)/SO{4) x SO{2). Left. Values of the DM relic density 
as a function of tudm for 2.0 TeV < F < 2.4 TeV fixing A = ^4 = de = A' = 0 and m^o = mH± = 
inriDM- The gray points correspond to the non-composite case with A = 0. We also show by dashed 
lines the 2cr CMB bounds on Right. <JoMp as a function of moM for points in agreement 

with the CMB bounds for F = 2.5 TeV and F —>■ 00 leaving all quartic and derivative couplings as 
free parameters. 


symmetry groups Q, it is possible that in the low energy effective theory one ends up 
with a triplet DM or even higher SU{2)l representations. As it is well known from non¬ 
composite models [42, 63], such cases are highly constrained and tend to require multi-TeV 
DM masses. In this section we will discuss how these features are affected by the inclusion of 
the dimension-6 operators induced by the strong sector. As discussed in section 3, we focus 
here on real representations, which correspond to ^ and simplify considerably the 

analysis. Also, after EWSB, the lightest component of the DM candidate from a complex 
representation usually is a charged field and do not provide a viable DM candidate [42]. 

The effective DM Lagrangian for any real DM multiplet is given by: 

£(2) = _ ^||^|2 _ A3 (^1 + |C|2|i7|2 

(3.24) 

-^\^mytQLHHR +ybQLHbn+ h.c.) , 

where the coefficients adi,d4, and A3 4 5 are taken to be real 0(1) parameters, the cou¬ 
plings Aj with i = 3,4,5 are allowed to vary in the window [—dvr, dvr] and /rg is the DM 
bare mass. The Pauli matrices are again denoted by = {cri, <72, <73} while the T matri¬ 
ces are the generators of the representation n of SU(2)i satisfying [rj,rj] = ieijkVk (see 
appendix A.3 for more details). 
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3.3.1 Generic composite triplet DM 

In order to provide some insight on the general behavior of a composite DM multiplet 
of dimension n > 2 of SU(2)£, we first focus on the case of the triplet. Later we will 
discuss how these results generalize to larger representations. In the triplet case, we use 
the following generators (in the spherical basis): 




r3 = 



(3.25) 


and ^ = {T~^,T^,T )^. Using eq. (3.24) with the above representation for the L matrices, 
we obtain the following scalar mass spectrum: 


'^‘dM — 

m‘^± 


2 2 A X 1 \ ^ 

t4 + ( 2 ^^''' 2 ■ 


(3.26) 

(3.27) 


Notice that in the triplet case, once the values of rriDM, A 4 and A 5 are given, the value of 
mrp± is fixed. 

In the non-composite case {F —)• 00 ), the triplet is exactly degenerate (except for small 
EW loop corrections, see [63]). In the latter case, either LEP constraints excluding the 
existence of low mass new charged particles or too strong annihilation into gauge (virtual) 
bosons^^ prevent the existence of viable DM candidate in the low mass range (i.e. for 
< w-h)- In file composite scenario, however, this is no longer the case since the 
dimension -6 operators induce a mass splitting, given by: 

- mjjM = ^ (^4 - As) . (3.28) 

Therefore, assuming |A 4 _ 5 | < 1, we can have itidm < 60 GeV imposing mrp± > 70 GeV, 
if E < 1.5 TeV. Larger values of F, however, will not generate sufficient mass splitting to 
allow for viable solutions in the low mass region. We will analyze this in more details in 
section 3.3.2. 

The effective Higgs-DM coupling in the triplet case is simply: 

=2~ ’ (3.29) 

/ v‘^\ v‘^ 

where A = As f 1Ag^ j -bdAs^ . (3.30) 


Unlike the doublet case, the DM-gauge boson couplings are only affected by one single 
coefficient ( 0 ^ 4 ) in eq. (3.24). While the contributions to the quartic vertices (^^V)(U^) are 
usually negligible: 

triplet 3 2g^ (^1 - adA^'^ T^T^W+W- , (3.31) 

the triplet case VV* is already dominant for moM ^ 50 GeV. 
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as they are always suppressed by j the triple vertex couplings can be enhanced 

at large moM- In addition let us emphasize that, for the triplet and higher representations, 
the coupling to the gauge bosons is always larger than in the doublet case because it 
scales like {v? — 1)^. In the case of the doublet, we saw that once moM — fnw, the 
(co-)annihilation into gauge bosons generate too large resulting in too low relic 

densities until ttidm ^ 500 GeV. This feature is only enhanced in the triplet case, since its 
couplings to gauge bosons are larger, as it is well known in the non-composite case [42, 63]. 
We also point out that, due to CP-invariance, Z decays to DM pairs are always forbidden, 
so the Z-width constraints do not apply to the case of real triplets. 

3.3.2 Scan results and experimental constraints 

In order to fully cover the allowed parameter space in a model independent approach, 
we once again scan over the effective Lagrangian couplings within the following ranges: 

20 GeV < rriDM < F , 

10-6 < Ah6 < 1, 

—dvr < A3 < dvr, 

— 1 < A 3 , A 4 , As < 1 , 

— 1 < C4, 02//, ^4 < 1 , 

—1 < ad4 < 1 • 


Again, we scan logarithmically over A’s ranges and linearly on all the other parameters. 




Xenon!T 

(projected) 

• LUX 

• BR(h -> inv.) 

• PAMELA 


120 


m^M(GeV) 


Figure 11. Generic composite triplet DM in the low mass range. Left Values of the DM-nucleon 
scattering cross-section as a function of tudm satisfying to 0.0941 < < 0.127. The channels 

contributing with the largest branching ratio to the annihilation cross-section at the time of freeze- 
out are depicted with different colors. Right: Experimental constraints on the low mass region. See 
text for details. 

The results of these scans for the low mass region are shown in the LH side of hg- 
ure 11, where we show the mass splitting and the DM-nucleus scattering cross-section as 
a function of rriDM- We take F = 0.8 TeV, since higher values of F will only reduce 
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the allowed region, due to the suppression of the mass splitting between the neutral and 
the charged component. We also show by blue (red) points the region where bb 

I'j'Oj’O W^*')W) dominates the annihilation cross-section in the early universe. We see 
that, for rriDM — 50 GeV, annihilations to (off-shell) gauge bosons become dominant and, 
for rriDM ^ 60 GeV, (cTann-'w) becomes too large, resulting in relic abundances below the 
GMB bounds. In the RH side of figure 11 we show the experimental constraints on the low 
mass region. As we can see, most of the parameter space is excluded by experimental con¬ 
straints, or theoretically disfavored for extremely low values of A (< 10“^). Only a small 
number of points within Higgs-resonance and W threshold and with A > 10“^ are still 
viable and will eventually be tested by XenonlT experiment. We conclude that, despite 
allowing for large mass splittings, the composite triplet scenario, alike the non-composite 
case, remains mostly excluded in the low mass region. 

In a way analogous to the doublet scenario, for moM S> mw, the annihilation cross- 
section into gauge bosons becomes suppressed enough to give rise to the right relic abun¬ 
dance. In this high mass region, it takes the following form: 


{(rann-v)\ToTo^vv - XX- 2 — (2|-Adr|^ + \Ml\^) , where 

olTTlTij^j^ 

A -2 


Mt ^ VtVt) ^ 2gl 


1 “d4 ^2 + 3 ^2 


^dl 


2^2 


f^DM 


— ( 1 “ + 


V 


7?2 o ™2 


A, 




Ml (TV ^ VlVl) ^ , 


(3.32) 


where we have taken rriDM — innrp+, since the small mass splitting can be neglected in the 
high mass region. Notice that in the large mass regime, we typically have (2m|)^) <C 

1, and for the F ^ v considered here we also have a^AV^ < 1. These two contributions 
to AIt will thus typically be small. 

In the LH side of figure 12 we show the high mass region results for F = 2.5 TeV and 
the non-composite case {F —>■ oo). As we can see, the lowest tudm values which generate 
the correct relic abundance are tudm ^ 1-7 TeV. Hence there are no possible solutions 
for F = 0.8 TeV, since our effective Lagrangian assumes ruDM < F. We also show the 
LUX and projected XenonlT constraints. As we can see, the large mass region has not yet 
been probed by LUX and will only be partially tested by XenonlT. From this plot, we see 
that while in the non-composite case a non-observation of signal in XenonlT would imply 
moM ^2.1 TeV, in the composite scenario this upper bound can be easily avoided due to 
the presence of the dimension-6 derivative operator. 

Let us emphasize that for such large triplet masses the Sommerfeld effect can be non- 
negligible [80]. A full treatment of the Sommerfeld enhancement within composite dark 
matter scenarios is however beyond the scope of this work. Let us mention though that, 
already in the non-composite case, one expects a non-perturbative enhancement factor of 
the effective coupling driving the annihilation cross section which is approximately constant 
away from the resonances and is about 1.6 in the triplet case [42, 80]. This implies an even 


- 28 - 










moM (GeV) (Q^y) 


Figure 12. Generic composite triplet DM in the large mass range: Values of the DM-nucleon 
scattering cross-section as a function of moM satisfying to 0.0941 < floMh'^ < 0.127. Right: Values 
of \adi\ in the aoMp — moM plane for F = 2.5 TeV. 

stronger lower bounds on muM in the high mass region. For instance, if the Sommerfeld 
enhancement [81] is included in the non-composite limit, one would obtain muM ^ 2.1 
TeV, instead of 1.7 TeV. 

The effective Lagrangian in the non-composite limit contains only two parameters (A 3 
and ttidm) and, as in the singlet case, the DM-nucleus cross-section value is fixed (for 
a given uidm) once the constraint on the relic abundance is imposed. The band seen 
in the LH side of figure 12 for F —)• 00 is simply due to the fact that we allow 
to vary within the interval [0.0941,0.127]. For the composite case, however, receives 
contributions from the dimension -6 coefficients A 3 and a^i- In this case, we can have larger 
or smaller values of aoMp (when compared to the non-composite case), depending on the 
values of adi- In the RH side of figure 12 we show the values of la^il in the aoMp — rnoM 
plane. We can clearly see that values of la^i | > 0.3 correspond to the regions where the DM- 
nucleus cross-section is either suppressed or enhanced with respect to the non-composite 
values. It is also interesting to notice that, while the low mass region requires F < 1.2 TeV 
in order to generate enough mass splitting, the high mass region requires F > muM ^1-7 
TeV. Hence there are no values of F for which there are solutions in both the low mass 
and the high mass regions. Furthermore, since the low mass region is already excluded by 
LUX, the triplet case can only provide a DM candidate for F > 1.7 TeV. 

3.3.3 Higher Representations 

Since the effective Lagrangian in eq. (3.24) is valid for any real DM representation, 
we can readily extrapolate the general conclusions obtained for the triplet case to higher 
dimensional representations. Since only odd representations contain a electrically neutral 
component, the possible representations higher than the triplet are the quintuplet (n = 5), 
the septuplet (n = 7) [63]. For any n-odd representation a mass splitting between the 
neutral (DM) and charged components is always induced by the dimension -6 operators 
proportional to A 4,5 in composite scenarios. We expect the allowed solutions in the low 
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mass region to be similar to the triplet case shown in figure 11. Consequently, in the 
case of higher representations, most of the low mass region is expected to be excluded by 
experimental constraints or theoretically disfavored. 

For the high mass range, given the results obtained for the doublet and triplet DM, one 
can expect that the threshold mass is the same in composite and non-composite scenarios. 
In the latter case, one can check that the total effective annihilation cross-section in the 
early universe has a given dependence in n [42, 63]. In the pure gauge limit it scales as 
(n^ — l){n? — 3)/n for a multiplet of dimension n so that for n = 5(7) one would have 
~ 4.3 (7.5) TeV. Again, one expect that Sommerfeld corrections will push these 
thresholds to even higher masses, see [80] in the non-composite limit. 

If we consider compositeness as a solution to the little hierarchy problem, we must 
require F < 3 TeV in order to obtain a fine-tuning of the order of 1% or lower [39]. 
Therefore, since moM < F, we see that the triplet is the highest DM representation allowed. 
Although there is still a small range of low DM masses allowed for all representations (50 
GeV < niDM ^ 60 GeV), these solutions usually require extremely small values of A^/i, 
which are unlikely to be generated in realistic models, where one expects X^h ~ 0(0.01 — 1). 
Gonsequently we conclude that Dark Matter representations higher than the triplet are 
highly disfavored in composite DM models. 

4 Conclusions 

The nature of DM is one of the greatest conundrums of our time. In spite of the fact 
that DM constitutes 85% of the total matter in the Universe, it continues to evade direct 
experimental observation. A reason for that may be connected to the fact that DM may 
be not elementary. Having this in mind we investigate composite multiplets that can arise 
in composite Higgs models and study the conditions under which they can be a suitable 
DM candidate. 

We considered a class of composite models where the only composite states present in 
the low energy effective theory are the Higgs and the DM, the first being a bi-doublet and 
the latter a multiplet of S'0(4), both pNGBs of a spontaneous broken global symmetry of a 
new unknown strongly coupled sector. We constructed and parametrized the most general 
effective Lagrangian up to dimension-6 operators under these general assumptions. We 
then checked if the different DM 50(4) multiplet candidates could account for all the DM 
assuming that the freeze-out mechanism is driving the relic abundance. For different DM 
representations under SU{2)l (or 50(4)) we derived the main DM observables (relic abun¬ 
dance, spin-independent DM-nucleon scattering cross-section, annihilation cross-sections) 
and imposed experimental constraints from direct and indirect detection experiments as 
well as LHC constraints from invisible Higgs decays. 

In the generic singlet DM scenario (arbitrary cosets), the DM-nucleon scattering cross- 
section can be significantly suppressed and can evade all current experimental bounds for 
"niDM ^ 100 GeV (except for a small discontinuity near the Higgs resonance), if T > 
0.8 TeV. However, once a specific model is considered (with a given coset), it can be 
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severely constrained by the DM observables. As an example, we discussed the case Q/"H. = 
50(6)/50(5). 

The composite doublet DM models differ from the singlet one due to tree level DM- 
gauge boson couplings and co-annihilations between DM multiplet components. There are 
two viable mass regions, one below the Higgs mass and one above ~ 500 GeV. Most of the 
points in the low mass region are either excluded by data or theoretically disfavored. For 
inn DM ^ 500 GeV, the annihilation into gauge bosons get suppressed enough to account for 
the DM relic abundance, just as in the IDM scenario. In the case Q/% = S'0(6)/50(4) x 
SO{2) presented in section 3.2.3, the dimension-6 operators allow for solutions only if T > 2 
TeV. The DM mass is then constrained to the range 600(550) GeV < itidm ^ 0.8(1.2) TeV 
for F = 2.1(2.4) TeV. 

The composite triplet DM scenario is also mostly excluded by data or theoretical 
considerations in the low mass region, despite allowing for larger mass splittings between 
neutral and charged components than in the non-composite scenarios. The large mass 
region will only be partially tested by XenonlT. We also mention that while in the non¬ 
composite triplet scenario a non-observation of a signal in XenonlT would provide an upper 
bound on moM, in the composite case this upper bound could be easily avoided. Finally, 
we have also examined higher representations and concluded they are highly disfavored in 
composite Higgs models if one requires small fine-tuning in the EW sector (or F < 3 TeV). 

We have shown that DM in various representations in composite models can reproduce 
the correct relic abundance and still be compatible with limits on the Higgs invisible width 
and from the non-observation of DM in direct and indirect detection experiments. These 
models will be further put to the test by future LHG, direct and indirect detection data, 
either excluding them completely or revealing some exciting new physics. 
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A Dimension Six Operators for Composite DM 


In order to compute the DM observables relevant for our results in section 3, it is 
essential to identify the relevant operators induced by the new strong sector. In sections 3.1, 
3.2 and 3.3 we presented the effective Lagrangians considered in our analyses. Here we 
will derive them in detail and justify (when necessary) our choice of operators. First we 
consider the SM dimension-6 operators. In addition to the ones listed in (eq. 2.4) there 
are other SM dimension-6 operators that are in agreement with the assumptions discussed 
in section 2 but were not considered in our analysis. Assuming minimally-coupled theories 
and a specific basis, the following dimension-6 operators can be induced at tree-level by 
integrating out heavy states with spin ^ 1 [48-53]: 


CHQ\ ' 
F 2 
CHui 

aw ig 
M 2 

a2W 9^ 

Mjgl 


{UR-i^UR) , ^ {dRj^dR) , 


(A.l) 


CHudi f TTCf^ 
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(A.2) 






(A.3) 


{D^w^,y{DpWf^y, 


,/2 


02 B g‘ 




Q 2 G ffs 

Mjgj 




(A.4) 


where D — {D^Uy a^H, Mp ~ gpF, and gp is the typical coupling of 

the Higgs doublet and SM fermions to the heavy resonances {gp < dvr). The coefficients at 
and Ci are 0(1) numbers. The operators in (A.3, A.4) have at least a suppression of 1/g'p 
compared with the operators in Additionally, the operator proportional to aw (A.3) 

contributes at tree-level to the S parameter, which is very constrained by the electroweak 
precision data [40]. For these reasons, we do not consider these operators in our effective 
Lagrangians. On the other hand, (A.l, A.2) have the same suppression factor (1/1^2) as the 
operators in Tg and should, in principle, be considered. However, these operators modify 
the gauge couplings to fermions which are in great agreement with the SM [40]. We can 
try to avoid this constraint by absorbing the universal part of the vertex corrections to a 
redefinition of the electroweak gauge fields [82, 83]. This restores the gauge couplings to 
their SM values, but it generates contributions to the electroweak precision parameters, 
which also are severely constrained by the experimental data. Therefore, all the operators 
in (A.l, A.2) are extremely constrained and we disregard them. We also neglect the 4- 
fermion operators since they do not affect the DM observables discussed here. Finally, 
operators that can only be generated at one-loop level, such as a^^yRH’^Wpv, are 
suppressed with respect to the operators in Tg and can also be neglected. 

In the remaining sections we compute the ^-dependent part of the effective DM La- 
grangian (T^^^) for specific DM representations. In order to be consistent and also to reduce 

^^Note that the suppression of the operators in (A.3, A.4) compared with Ce is larger as the theory is 
close to the strong-coupling limit, i.e. when the typical coupling is Qp ~ 47r. 
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the number of free parameters, we assume in the DM sector the same simplifications pre¬ 
viously required in the Higgs sector, i.e. 


• we neglect the suppressed operators; 

• from the electroweak precision data [40] we know that the c’s in (A.l, A. 2 ) must be 

very small <C 1). Here we assume that the analogous coefficients in the DM 

sector are also suppressed. 

A.l Singlet DM 

The case of singlet DM is trivial since ^ can only couple to operators which are singlet 
under SU{2)l x U{1)y- The only operators (up to dimension 6 ) containing two DM fields 
and the Higgs field contributing to the scalar potential are given by 

(A.5) 

The dimension -6 derivative interactions, corresponding to the effective CCWZ, must pre¬ 
serve the symmetry 50(4) C Ti, hence the only allowed operators are 

. (A. 6 ) 

The last two operators can be eliminated after the following field redefinitions 

+ (A.7) 

then, 

(Dpif)t O'-if + ^ D^H + 

From the above result we see that choosing the coefficients a and a' properly, the opera¬ 
tors and (D^H)^ D^H can always be eliminated. Therefore, we obtain the 

following effective Lagrangian for the singlet DM field 

i:® = - \pie - y (i + 

- (ytOLii^fR + VbQL^bf- + h.c.) , (A.8) 

where ^4 is real as we assume a OP-even Higgs. 

A.2 Doublet Dark Matter 

The doublet DM, ^ ~ (2,2), is quite involved since in this case ^ has a non-zero 
hypercharge. However, the number of possible operators can be greatly simplified imposing 
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50(4) invariance and using suitable field redefinitions. In order to simplify the notation, 
we use the following 50(4) bi-doublets to construct the invariant operators 

^ = and = (f ,0 • (A.9) 

The above 2x2 fields transform under 50(4) as 

and , (A. 10) 

so Tr [<I>'l’<b] is a 50(4) singlet. 

For the scalar potential we can construct the following SU(2)l x U{1)y invariant terms 
involving 2 powers of ^ 

dimension 2: Tr ; (A.11) 

dimension 4: Tr Tr , Tr , 

r 1 2 r 1 r 1 

Tr 4>'^<l)^fT^ , Tr Tr <l)'l’4>go-^ ; (A.12) 

r i^r+i r 1 r 

dimension 6: Tr Tr , Tr Tr , 

rir i 2 rir ir 1 

Tr <b'^<b Tr <l)^<bg(T^ , Tr 4)'I’<I> Tr 4>'I’<1)^ Tr <I>^’<b^cr^ . (A.13) 

Notice that all the operators involving explicitly break 51/(2)^. However, they are 
allowed in the scalar potential since it is generated by the Yukawa and gauge couplings 
between the SM and the strong sector, which violate the 50(4) symmetry. All the other 
possible operators (including the triplet-triplet contractions) can be reduced to the ones 
above. Writing the above operators in terms of ^ and H we obtain the following scalar 
potential 

y{H, 0 = mIICP + As (l + + A4 (l + 

+ y + . (A.14) 

In the doublet case there are two dimension-6 Yukawa-type operators involving the 
singlet-singlet and triplet-triplet operators 

Wuk(<&,0 = {ytQLHHR + yi,QLHhR + h.c.) 

+ {ytQLSHHR - ybQi^HbR + h.c.) 

+ ^ + yti'^ai^^QiSHtR + h.c.^ . (A.15) 

In order to avoid too many irrelevant parameters in the potential we take dg = which 
we do not expect to affect any of our results. 
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Finally, the dimension-6 50(4) invariant derivatives operators, corresponding to the 
effective CCWZ Lagrangian, can be constructed from the operators listed in (A.11, A.12, 
A.13) 


Tr 

Tr 

Tr 




Tr 


, Tr 


Tr 


, Tr 

, 

Tr 







1 2 

, TV 


2 

, TV 


Tr 


, Tr 


Tr 



- Tr 




(A.16) 


The first three operators can be rewritten in terms of the others if we consider suitable 
field redefinitions 


H ^ H + aTi 
i +a'T t 




H + bTi 




^ + 6 Tr 




I, 

H, 


(A.17) 


where a and a' are chosen to cancel the first two operators and b is chosen to cancel the 
third one. Furthermore, we notice that 


Tr 




1 2 


-Tr 




1 2 


= d^Tr 




Tr 




-Tr 




(A.18) 


The above operator can also be eliminated through the field redefinition 


H ^ H + cTr 




e, e^^-cTv 




H. 


(A.19) 


with the coefficient c suitably chosen. Therefore, after the above field transformations we 

obtain only 4 independent operators, resulting in the following effective CCWZ Lagrangian 
15 


c. 


( 2 ) 

CCWZ 


= D^^H + 


adi 


+ ^TV 

^ J?2 






Tr 




+ 


2 F 2 

«(i3 


Tr 




Tr 


Tr 




+ 


adA 

i ?2 


Tr 




-Tr 


i ?2 


+ Tr 




2 

(A.20) 


Finally, rewriting the operators in terms of the H and ^ fields and combining (A. 14), 
(A.20) we have replaced the operator 

Tr " -h Tr " by (Tr [4>tD^4-e] -t Tr [4>JL»^4>] ) " , 

since they only differ by the operator proportional to ad2- 
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(A. 15) and (A.20), we obtain the effective Lagrangian for the doublet DM field 


- As (l + - A4 (l + 

hY + h.c. 


2 


- — ( l + ^\H\'^ 


i? 2 l 


d^\H\^d>^\C\^ 


(A.21) 


+ ^(FtD^^ + h.c.) 


+ 


Qd 2 

i ?2 

F 2 L' 

di 


I 1 \ I ®(i3 

+ h.C.) 


dfj, + h.c.^ 




- [ytQiHHR + VbQLHbR + h.c.) 


i ?2 

d% 

L' 


{vtQLdHHn - ybQLdHbn) + yb^^'^d^QiaH^R + yti"^a^'^QiaHtR + h.c. 


where and dg are real as we assume that the Higgs doublet is CP-even. We have verified 
that the imaginary parts of A 5 and A 5 have almost no impact on the allowed parameter 
space. Therefore, for simplicity, we take these complex parameters as real. 


A.3 Real Representations 

Here we discuss the case where ^ is a real multiple!, which corresponds to the 50(4) 
representations (n, 1). The possible operators involving a real ^ multiple! are greatly re¬ 
duced due to the identities 

n + h.c. = rp'^e = o. (A.22) 


Furthermore, since ^ is a singlet under SU{2)r, the SO{4:) invariance implies that ^ can 
only couple to SU{2)r singlets. The only possible operators contributing to the scalar 
potential with two powers of ^ are 


dimension 2: Tr 


<[.*3-4 i 


(A.23) 


dimension 4: Tr 




(A.24) 


dimension 6 : Tr 
Tr 




Tr 




,T4 




Tr 








(A.25) 


where = cr^ ± ia"^. The dimension -6 derivative operators invariant under ^©(d), which 
parametrize the effective CCWZ Lagrangian, are given by 


Tr 


> l^pTr [p^4>1'P/,$ 


Ti- 


4>tP/,4>l and Tr [4>VP/,4>1 


(A.26) 
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As before, the first two operators can be eliminated through the field rescalings: 

H ^ H+ a\C\^H andC^C +a (A.27) 

Therefore, writing the above operators in terms of the H and ^ fields we obtain 

- As (l + 

- {r\F} iH^a^HH^a^H - |r\F} H (A.28) 

- [ytQLHHn + ybQiHbR + h.c.), 

where once again we take (i 4 to be real as we assume a CP-even Higgs. 
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